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We study the non-radial oscillations of relativistic neutron stars, in particular the (fun- 
damental) f-modes, which are believed to be the most relevant for the gravitational 
wave emission of perturbed isolated stars. The expected frequencies of the f-modes are 
compared to the sensitivity range of Mario Schenberg, the Brazilian gravitational wave 
spherical detector. 
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1. Introduction 

Why study oscillations of neutron stars? When it comes to sources of gravitational 
waves, binary systems of compact objects are among the most promising sources. 
But after a binary system inspiralls and merges, we are left with one perturbed 
compact object, which will somehow radiate away its extra energy to achieve a new 
equilibrium state, including gravitational radiation. And this final compact object 
could even be a super-massive neutron star, depending on the initial components 
of the binary system and the equation of state of the neutron star. Also single 
(isolated) neutron stars can be perturbed due to the interplay between the fluid 
of which the star is composed, its solid crust and magnetic field, for instance. Of 
course, this is a very complicated problem, beyond the scope of our present paper. 
But any generic perturbation with I > 2 can lead to the generation of gravitational 
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waves. The stellar modes of oscillation are classified according to the restoring force 
acting on themUJ They form different families of modes, with separated frequencies. 
For a non-rotating star (without magnetic fields) the usual classification gives: 

• /-modes (fundamental) : frequency is proportional to the square root of the 
mean density of the star. These are favoured for the gravitational wave 
emission. 

• p- modes (pressure): pressure is the restoring force, and the frequencies are 
higher than those of the f-modes. 

• g- modes (gravity): buoyancy is the restoring force, frequencies are lower 
than those of the f-modcs. 

For the spherical antennas, such as the Brazilian Mario Schenberg and the Dutch 
Mini-GRAIL ®, the /-modes would be particularly interesting due to the possibility 
of detection inside the 2.8-3.4 kHz frequency range, where they have their sensitivity 
bands. 

This paper is organized as follows: in section[2]we present our equilibrium stellar 
model and in section [3] we review the perturbation equations for the /-modes. We 
discuss the numerical method and present our results in section |4] and conclude with 
our final remarks in section [5] 

2. Equilibrium Star 

We begin with the spherically symmetric metric 

dS 2 = -e v ^dt 2 + e x ^dr 2 + r 2 (d6 2 + sin 2 6d<f> 2 ), (1) 

and the perfect fluid stress-energy tensor 

T"" = (e +pg> tu , (2) 

where e is the total energy density, p is the pressure and m m is the fluid 4-velocity 
to write the TOV equations 

m = Airer , v = — — — , p = -v — - — . 3) 

r (r — 2m) 2 

Finally, we use a polytropic equation of state p = Kp T , where p is the rest-mass 
energ density of the fluid ^1. 

3. Polar Perturbations 

The linearized Einstein equations 

S (R„u - ^uRj - 87r<5T M „ and 5{T^) = , (4) 

with the following Ansatz (polar perturbations) for the perturbed metric tensor 
dS 2 = -e v {l + /HoYie^dt 2 - 2iujr e+1 H 1 Y^e lu ' t dtdr + 

+ e A (l - /H Yle iut )dr 2 + r 2 (l - r l KY^e lU3t )(dQ 2 + sin 2 9dtf) , (5) 
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and for the perturbation of the fluid in the star SI (Lagrangian displacement) 



lead to a 4th order system: 
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where X = u> 2 (e + p)e~ v / 2 V - ^. e {u ' x)/2 W + |(e + p)e v l 2 H and i? is given 
by another algebraic relation EE1 



4. Numerical Method and results for the /-modes 

We look for solutions that are regular at the origin r — 0, have vanishing perturbed 
pressure at the surface r — R and describe only outgoing gravitational waves at 
infinity r — > oo. 

Our numerical code follows closely the prescriptions of Ref. |6] We use as an 
initial guess the Newtonian frequency of the f-modes: 




21(1 - 1) / M \ 
2^+1 \~Rz) ' 



and we solve the eigenvalue problem with a modified shooting method. By combining 
at R/2 the two linearly independent regular solutions at r = with the three linearly 
independent solutions with vanishing pressure (X(R) = 0) at the surface we create 
the initial conditions for solving the Zerilli equation outside of the star, imposing 
then only outgoing waves at infinity. 

In Figs. [1] and [5] we show the perturbation functions inside and outside the star 
for a typical case (with K = 100, T = 2 and p = 1.28 x 10~ 3 in the equation of 
state, given here in code units, i.e. c = G — M sun = 1). 

Typical values for the complex frequency uj = 27r/ + ijr of the /-modes for a 
neutron star have / = 1.5 — 3 kHz and r = 0.1 — 0.5 The values obtained for 
typical star shown in Figs. [T] and [2] are / = 1.5795 kHz and r = 0.2987 s. 

It has been observed that / and r can be fitt ed by simple expressions in terms 
of the total mass M and the radius R of the star p or a sequence of stars with 
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our typical polytropic equation of state and increasing central density, we produced 
the following fit: 



9 [W 1 M 3 

f = 7.36 x 1CT 2 + 55.80a/ and - = — j- 
V it J r it 4 



9.91 x 10" 2 -0.33 ( — 



(8) 



which can be seen in Fig. [31 together with two other fits from Ref. [51 that used 
realistic equations of state and equations of state containing strange quarks ^. 
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Fig. 1. Behavior of the perturbation functions inside the star. 





Fig. 2. Perturbations outside the star, given in terms of the Zerilli function Z and its derivative, 
as functions of the tortoise coordinate r t . 
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5. Observational Consequences and Final Remarks 

Ideally, after observing the /-mode frequencies from neutron star emitting gravi- 
tational waves, we will be able to solve the initial value problem and determining 
the total mass and radius of the star. With this information, we will gain precious 
information on the neutron star equation of state, which is still largely unknown, 
and object of investigations from different branches of physics. 

From the observational point of view, the /-modes are the most relevant ones, 
expected to be favoured for the gravitational wave emission. Our fits for the fre- 
quencies, even though obtained used only a very simple polytropic equation of state, 
already present the expected behavior and typical values. The Mario Schenberg and 
Mini-GRAIL spherical antennas may observe these neutron star /-modes in a near 
future, if some of these modes fall inside the 2.8-3.4 kHz frequency range, where 
they have their sensitivity bands. 
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Fig. 3. Numerical fits for the frequency and damping time of the /-modes. 



